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Abstract. We study the effects of random scatterers on the ground state 
' of the one-dimensional Lieb-Liniger model of interacting bosons on the unit 

^ , interval in the Gross-Pitaevskii regime. We prove that Bose Einstein condensa- 

■ tion survives even a strong random potential with a high density of scatterers. 

The character of the wave function of the condensate, however, depends in an 
essential way on the interplay between randomness and the strength of the 
two-body interaction. For low density of scatterers or strong interactions the 
wave function extends over the whole interval. High density of scatterers and 
weak interaction, on the other hand, leads to localization of the wave function 
in a fragmented subset of the interval. 



1. Introduction 

While the effects of random potentials on single particle Schrodinger operators 
P] and ideal Bose gases 2-4 are rather well explored, the present understanding of 
such effects on many-body systems of interacting particles is much less complete. 
In recent years, however, many papers concerning the interplay of Bose-Einstein 
Condensation (BEC) and disorder have appeared of which references [S]-[T2] are 
but a sample. 

We present here results on a model that is in a sense the simplest one imaginable 
where this interplay can be studied by rigorous mathematical means. This is the 
one-dimensional Lieb-Liniger (LL) model [17 of bosons with contact interaction in a 
'flat' trap, augmented by an external random potential that is generated by Poisson 
distributed point scatterers of equal strength. We study the ground state and prove 
that, no matter the strength of the random potential, BEC is not destroyed by the 
random potential in the Gross-Pitaevskii (GP) limit where the particle number 
tends to infinity while the coupling parameter in a mean-field scaling stays fixed. 
The character of the wave function of the condensate, however, depends in an 
essential way on the relative size of the three parameters involved. These are the 
scaled coupling parameter 7 for the interaction among the particles, the density of 
the scatterers v, and the strength a of the scattering potential. All these parameters 
are assumed to be large in suitable units. 

Our main results an be summarized as follows: For a ^> 1 we identify three 
different "phases" of the condensate. For large 7 3> v 2 the condensate is extended 
over the whole trap (the unit interval in our model). A transition from a delocalized 
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to a localized state takes place when 7 is of the order i/ 2 , in the sense that for 7 <C v 2 
the density is essentially distributed among a fraction A <C 1 of the v 3> 1 intervals 
between the obstacles. For 7 3> f/(ln z^) 2 we still have Af ^> 1, but for 7 ~ ^/(ln r/) 2 
the fraction of intervals that are significant occupied shrinks to 0(y~ v ). We stress, 
however, that in all cases there is complete BEC into a single state in the limit 
when the particle number tends to infinity. 

Details of our analysis are given in the preprint larXiv :1207.7054. 

2. The Model 

The model we consider is the Lieb-Liniger model of bosons with contact inter- 
action on the unit interval but with an additional external random potential V u . 
The Hamiltonian on the Hilbert space L 2 ([0, l],dz)® eN is 

N 

(1) H = J2 {-dl + V u {Zi)) + ^ - z i) 

i=l i<j 

with 7 > and Dirichlet boundary conditions. The random potential is taken to 
be 

(2) 

3 

with a > independent of the random sample u> while the obstacles {z^} are 
Poisson distributed with density i/>1, i.e., their mean distance is v . 
There are several reasons for studying this model: 

• It is the simplest model of its kind. 

• BEC in the ground state can be proved in a suitable limit, for an arbitrary 
nonnegative external potential. 

• V u is simple enough to allow a rigorous analysis of the condensate. 

We remark also that the case 7 = 0, a — 00 corresponds to the Luttinger-Sy model 

0- 

Once BEC has been established the main question concerns the dependence of 
the properties of the condensate on the three parameters 7, a and v. 

3. BEC in the Ground State 

A basic fact about the model ([1} is that for fixed 7, a and configuration u) there is 
complete BEC in the ground state in the sense that the 1-particle density matrix/iV 
converges to a one dimensional projector as N — > 00. The corresponding wave 
function of the condensate is the L 2 -normalized minimizer of the Gross- Pitaevskii 
( GP) energy functional 

(3) £ GP M = f {W{z)\ 2 + V u (z)\i/,(z)\* + ( 7 /2)|V^)| 4 } dz. 

Jo 

We want, however, to consider is, a and 7 large. Hence it is important to estimate 
also the rate of the convergence of the 1-particle density matrix as JV — > 00, in 
dependence of these parameters and of the configurations u. 
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3.1. The Proof of BEC (sketch). The proof of BEC is based on energy bounds: 

• An upper bound to the many-body ground state energy Eq m by taking 
ipQ N as a trial function for H where tpo is the minimizer of the GP energy 
functional, normalized so that HV'olh = 1- This gives 

(4) £ QM < Ae 

where eo = £ GP [ipo] is the g.s.e. of the GP functional. 

• An operator lower bound for the many-body Hamiltonian H, up to con- 
trolled errors, in terms of the 1-particle mean-field Hamiltonian 

(5) h = -dl + V u (z) + !\Mz)\ 2 ~ (7/2)/|V>o| 4 

which has ?po as ground state with energy eo. 

BEC follows from the upper and lower bounds and the fact that there is an energy 
gap between eo and the next lowest eigenvalue, e%, of the mean-field Hamiltonian 
h: 

Let 

(6) A = (* |ot(V>o)o(^o)|*o> 

be the occupation number of the GP ground state ipo in the many-body ground 
state "to- Then the energy bounds give 

(7) A e + (N - N )e 1 - o(l)Ne < E$ M < Ne 

where the o(l) factor depends only on 7 and N. This implies an upper bound for 
the depletion: 

(8) (l-^)<o(l) 60 



N J ei — eo 

More precisely, the estimate on the depletion of the condensate is 
THEOREM 1 (BEC). 

(9) (i-^2.) < C^^N- 1 ^ miniy/ 2 , 7}. 



N J ei - e 

We must, however, also consider the dependence of the energy gap, ei — eo, on 
the random potential and the parameters. 



3.2. Remark. If N<k denotes the occupation of the k lowest eigenvalues of §5§ 
with energies up to e&, then we have 

(io) (1 ^) < c^r'/^h 1 ^}. 

V A J e k -e 

For finite A this estimate may be more useful than (|9|) because even if k N, 
efc — eo can be substantially larger than ei — eo. Thus the right side of (flO|) can be 
< 1 even in cases when ((9|) contains no information because the right-hand side is 
> 1. 



4 



ROBERT SEIRINGER, JAKOB YNGVASON, AND VALENTIN A. ZAGREBNOV 



3.3. The energy gap. Consider a one-dimensional Schrodinger operator —d 2 + 
W(z) on the unit interval with a nonnegative potential W and Dirichlet boundary 
conditions. 

LEMMA 1 (Gap). Define n > by n 2 = n 2 + 3 /„* W(z)dz. Then 

(11) ex - e > rj In (l + Tre" 2 '') 

The proof is based on a modification of a result of Kirsch and Simon [19] . 
that involves the sup norm of W instead of the integral. In our case rj = r\ u = 
\J 7T 2 3m w <7 + 37 where to w is the number of obstacles in [0,1], that is almost 
surely equal to v in the limit v — > 00. 

For large am^ the bound is certainly not optimal, in fact, in this case one expects 
e x - e ~ (crm^y 1 . 

3.4. The Poisson distribution of the obstacles. Let z" < z% < ■ ■ • < 

denote those random points which lie in the unit interval [0, 1]. The lengths £i — 
zf +l — zf are independent random variables with distribution 

(12) dP v {l) = ve- iu dl 

and we are considering the case v ^> 1. 

The average length of an interval free of obstacles is v~ l and with probability 
one, m^/v —y 1 and y^^ij = 1 for v 3> 1. Combined with the estimate on the 
energy gap this implies that the depletion of the BEC is uniform in the L p norm 
on sample space for any p < 00. 

3.5. Ideal vs interacting gas. While the average length of an interval is v~ Y there 
is, with probability one, a unique largest interval of length ~ v~ x \vlV. If there is 
no interaction, i.e., 7 = 0, and a —> 00 (Luttinger-Sy model), then the ground state 
wave function will be localized in the largest interval for kinetic energy reasons, 
with energy ~ ^ 2 /(lni/) 2 . 

The question is how the character of the wave function of the condensate, ipo, 
changes when the interaction, i.e., the term (7/iV) J \ipo(z)\ 4 dz comes into play. 



4. The Mass Distribution of the Condensate 

4.1. A limit theorem for the GP energy. Our first result is that the energy 
becomes deterministic in an appropriate limit. 

Let e w (7, cr, v) denote the GP energy with the random potential and eo(7, v) 
the energy for a — 00, averaged over to. 

THEOREM 2 (Convergence of the energy). Assume that v — > 00, a — > 00 and 
7 — y 00 in such a way that 

(13) 7 > — and a > 



(ln^) 2 1 + In (1 + ^/7) 

Then, for almost every sample to, 

(14) lta^ = l. 

e (7,^) 



DISORDERED BOSE EINSTEIN CONDENSATES 



WITH INTERACTION 



5 



4.2. Comments on the proof. The proof has two interrelated parts: 

• Comparison of e u (j, a, v) with e w (7, oo, v). 

• Comparison of e w (7, oo, v) with the deterministic eo(7, v). 
Calculations are conveniently done in a grand canonical ensemble, introducing a 

chemical potential \i that determines the optimal repartition 

(15) nw«(^rv 2 -^] + 

of the condensate mass in intervals of length I between the obstacles. The lengths 
are distributed according to dP v {€) and the normalization requires 

(16) v J n{i)dP v {£) = 1. 

4.3. Discussion of the GP wave function for large a. For large a the average 
number of intervals with non-zero occupation numbers is given by 

/•oo 

(17) v I dP v (t) = e-™'^v. 

Since v is the total (average) number of available intervals, 

(18) A := e'™!^ < 1 

defines the fraction of them which are occupied. The normalization requires 

(19) 1 ~ ^e-WVF 

7 

so A is determined by 7 and v via relation 

A 7 



(20) 



(In A- 1 ) 2 v 2 



We can now distinguish the following limiting cases: 

• If 7 > v 2 then by (gDJ) we get A — > 1, i.e., all the intervals are occupied 
(derealization). The chemical potential satisfies [i ~ 7 in this regime. 

• If 7 ~ v 1 then A ~ 1, but A is strictly less than 1 (transition to localization). 
Again we have fi ~ 7. 

• If 7 <C v 2 then A <C 1, i.e., only a small fraction of the intervals are occupied 
(localization). The relation p^|) implies fi ~ ? for the chemical potential. 

• If 7 ~ ^/(lni/) 2 then by ([2H)l the fraction A becomes 0(l/v), i.e., only 
finitely many intervals are occupied. In this latter case, \i ~ 71^ ~ z/ 2 /(ln j/) 2 , 
which corresponds exactly to the inverse of the square of the size of the 
largest interval. 

In particular, Xv 3> 1 only if 7 3> ^/(ln^) 2 , and hence this condition guarantees 
that many intervals are occupied. In this case the law of large numbers applies and 
hence the energy becomes deterministic in the limit. 

If \v = O(l), on the other hand, the value of ecj(7, a, v) is random. This shows, 
in particular, that our condition on 7 is optimal. 

Also the condition a 3> v/{l + ln(l + i> 2 /j)) can be expected to be optimal. It 
can be rephrased as la 3> 1, where £ is the (weighted) average length of occupied 
intervals 

/>oo 

(21) l = v / dP v (£)£n(£) 

Jo 
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with n(£) the optimal repartition of the mass. A simple calculation shows that 
^~i/- 1 (l + In(l + i^/7)). 

5. Conclusions 

• BEC in the ground state of the interacting gas in the GP regime can survive 
even in a strong random potential. The character of the wave function of 
the condensate, however, is strongly affected by the interaction. 

• A random potential may lead to localization of the wave function of the 
condensate in subintervals. The interparticle interaction counteracts this 
effect, however, and can lead to complete derealization (the condensate 
extends over the whole unit interval) if the interaction is strong enough. 

• In terms of the interaction strength, 7, and density of scatterers, v, the 
transition between localization and derealization occurs in the model con- 
sidered when 7 ~ v 2 . For 7 < v / ih\v) 2 a "third phase" occurs where the 
condensate is localized in a small number of subintervals. 
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